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Abstract 

In previous work the parameter of glassiness was introduced to distinguish between a liquid and 
a glass, using a formal analogy with the quantum Bose system. The glassiness is defined in such 
a way that it is unity in a frozen system and less than one in a liquid. In the present letter we 
revise first the results obtained for the glassiness in a hard sphere liquid as a function of the 
density. Then we investigate the influence of an attractive potential by obtaining the glassiness 
as a function of the density, temperature and the attractive tail when a square well potential is 
added to the hard core. 
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The transition from the liquid state to the glassy state can not be seen as a change in an 
equilibrium parameter, as when dealing with the liquid-solid phase transition. Transport 
parameters such as the diffusion constant of an external probe and the viscosity are the 
main criteria for distinguishing the glassy state from a simple classical liquid j 12 1 3 j . Some 
years ago an analogy between a simple liquid and a Bose quantum system was used to 
suggest an equilibrium parameter, which distinguishes between a simple liquid and a glass 
[1] . That parameter is based on the difference in the symmetry under permutation of the 
classical liquid and the glass. In both states the system is invariant under permutations 
but still there is an important difference in the physical effectiveness of the symmetry. In 
the liquid state configurations which differ by permutations are connected by trajectories 
in configuration space, which do not have to cross high energy barriers. In the glass, on 
the other hand, such trajectories have to cross extremely high barriers. In quantum 
systems the Bose condensed fraction is the parameter that measures the effectiveness of 
the symmetry under permutations. In the following we present the analogy between a 
classical system at equilibrium and a quantum Bose system and show that a condensed 
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fraction can be denned, which is non zero in the liquid state and zero in the glass, 
where the symmetry under permutations is not relevant. The glassiness , T, is defined as 

r = i-e (i) 

We begin with a set of N coupled Langevin equations for a classical system of volume 

V. 

Where W is the interparticle potential and T)i is a noise experienced by the j'th particle, 
obeying 

(>fc)=0 and ( % (t)» 7i (t')>=X>M(*-* / )- ( 3 ) 

The above provides a framework which allows a natural analogy between the equilibrium 
state of the classical system and the ground state of a corresponding quantum system. 
The Fokker-Planck equation for the distribution ,V, that corresponds to the above set 
of Langevin equations has an equilibrium solution V s = exp[— W/HbT], where ksT = 

1/2 

-D/27 [3]. The standard transformation V = V S ^, yields an imaginary time Schrodingcr 
equation for ^, 

St = < 4 > 

The Hamiltonian above is a non-negative definite Hcrmitian operator and its only state 

1 /2 

with eigenvalue zero is the ground state, H>g = V s ■ Because it is symmetric under 

1/2 

permutations it is also the ground state of the Bosonic reduction of TC. The state V s 
is thus a natural Bose "quantum" ground state that corresponds to the classical system 
and the question of the value of the condensed fraction becomes a legitimate question for 
the classical system, which as discussed above has bearing on the question whether the 
system is a liquid or a glass. The purpose of the present article is to study the glassiness 
in the liquid state for the hard sphere system and then consider the effect of adding an 
attractive interaction. For a potential energy which is the sum of two body interactions 
4>{ r i ~ r j)i the ground state wave function is given by 



$G = ex P 



N 

- ^(Xi-Xj) 

ijij=l 



(5) 



where f3 — l/fc^T. Some basic algebra was used in ref. [1] to obtain the condensed 
fraction, 

_ Q N+ i GQ3/2) 
? VQ N G(/?) ' 1 ' 

where Qn is the configurational partition function of a system of N particles and the 
function G(a) is defined by 



G(a) — < exp 



-2a J dx>(x') j o(x'; 



(7) 



where p(x) — Y^iLi ^( x — x i) an< i (" ' ) denotes thermal average with respect to V s - The 
volume limit is implied in eq.([7]) . For the hard sphere case the G factors cancel, so that 
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the condensed fraction can be easily expressed in terms of the chemical potential. We 
begin with the following relation 

^=«- W -> (8) 
where F ex is the configurational part of the free energy and is given by 

N 



fr — 



rfM (9 ) 

h \ Nk B r L J r,' ' (9j 



where rj is the packing fraction which is given in terms of the density p and the hard 
sphere diameter a as rj — ira 3 p/6 and P is the pressure. The excess chemical potential 
associated with this free energy is 

The glassiness of the hard sphere system, just like all other thermodynamic properties is 
independent of the temperature and is given by 

'Qx+A IV' 



iN 

so that 

Tffs = 1 - exp(-(3(i ex ) =l-exp 




(12) 



The equation of state for the hard sphere system has been obtained in the literature by nu- 
merous analytical approximations such as Percus-Yevick (PY) [5] , the Carnahan-Starling 
(CS) [6], the virial expansion (VEx)[7] etc. All of the above yield analytic expressions for 
the dependence of the pressure on the packing fraction which lead directly to analytic 
expressions for the glassiness. For the virial expansion the expression is trivial, of course. 
The Percus-Yevick and Carnahan-Starling produce respectively the following equations 
of state: 

(PV\ _ 1 + V + if 

y-Nkgrjry-— w { ' 

UwL = a-.) 3 ' (13b) 

which lead to the following expressions for the glassiness respectively, 

T PY =1-(1- V )eii 1 - 1 ^ 1 -^ (14) 

and 

_(i-3rf)r L 

T cs = 1 -e l^W . (15) 

Figure [T] presents the glassiness of the HS system, following from the approximations 
considered above. As can be seen from figure [1] it is not easy to distinguish among the 
various approximations for the glassiness. In figure we take the VEx glassiness as a ref- 
erence and present the relative departure of the PY and CS glassiness from the reference. 

We note that the results obtained above for the classical hard sphere system can also 
be used to estimate the condensed fraction in liquid H e 4 under various condition. This 
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Fig. 1. Classiness of the HS system in 3D. 




n 

Fig. 2. Relative differences of the glassincss calculated from the three different EOS. 

observation is based on the work of Penrose and Onsager [5] who used the classical hard 
sphere weight function as a model of the ground state function of liquid He 4 . They used 
that wave function to obtain the condensed fraction at 0.28 of the close packing density 
by employing the low density expansion to its lowest order. Our present calculation goes 
beyond the lowest order and supplies therefore, the condensed fraction in He 4 in a wider 
range of densities accessible under pressure, within the Penrose Onsager model. 
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The hard sphere system is a canonical model used in the study of simple liquids. Yet, 
realistic potentials have an attractive tail, which must affect the physical properties of 
the system including the glassiness. The effect of an attractive tail seems unclear. On 
one hand we expect it to counteract the repulsive potential and thus the introduction 
of an attractive tail seems to lead to a reduction of the glassiness. The addition of an 
attractive potential is expected on the other hand to increase the order in the system 
and is expected thus also to increase the glassiness. The overall effect seems to be density 
dependent. At low densities the first effect is expected to dominate, because introduction 
of an attractive tail will not increase the tendency to order and the glassiness will decrease. 
At high densities the latter effect will prevail. To get a basic idea of what is going on we 
add to the hard core a soft attractive potential. The two body potential is thus 

</> = 4>hs + 4>sp- (16) 

To first order in the soft potential the glassiness is given by 

r = T HS + (l-T HS )A S P, (17) 

where the dimensionless quantity, Agp, is the first correction to the condensed fraction 
divided by the condensed fraction of the hard sphere system, which is just the first 
order correction to In £. We can calculate the first order correction to In £ by considering 
separately In (Q,i+i/VQn) = —fipex and lnG(/3/2)/G(/3). The first order correction to 
the first term was obtained in the past [9110] and is given by 



Ac 



2 dp 



P 2 9 H s( r , p)(j>sp{r)dr 



(18) 



where g HS (r, p) is the radial distribution function of the hard sphere system at density 
p, which is assumed to be given. We obtain the first order correction to the second term 
to be given by 

A| P = /3p 2 [ 9 HS (r,p)<t>Sp(r)dr. (19) 



Thus, to first order in the soft potential the total correction is given by 

Asp = — J -Q-p 0sp(r)dr. (20) 

A common idea used in the theory of classical liquids which is considered to give a fairly 
accurate description of the liquid is to use the simplest attractive potential which is the 
square well potential defined by 

' 0, |r| < a 

"MrU-i = <(-e, <t< |r| < \<j (21) 
0, Act < |r|. 

For the SW potential the expression PO)) can be simplified: 

-e(3p 2 [ x ° dg HS (r,p) 



-dr. (22) 
dp 
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In order to evaluate this correction, we use an expansion in the density of the pair 
distribution function, 

oo 

g HS (r, p) = 1 + J2 H n (r/a)(pa 3 r for r > a. (23) 

n=l 

This density expansion was obtain on the basis of some analytical work and MC inte- 
grations up to six's order [IT] and is hereby used to calculate the integral in (|22j) (There 
are other forms that give the radial distribution function and the structure factor as a 
function of the density |12ll3ll4ll5j but the form given by (|2"3"|) proved to be the most 
convenient for our purpose). This results in explicit dimensionless coefficients a n {\) for 
n = 1 — 6 in the expansion of A$w m the packing fraction 77, 

12 6 

71=1 

where T* — l/fle is a reduced temperature. The correction, Asw, for T* = 1 is pre- 
sented in figure [3] as a function of the packing fraction for different well widths and in 
figure 2] as a function of the well width for different packing fraction values. We see that 




n 

Fig. 3. Asw as a function of the packing fraction, T* = 1. 

the results presented in the figures agree with our previous expectations. The correction 
decreases first as the density is increased reflecting the fact that the attractive potential 
counteracts the hard core. Then the correction starts to increase as the density is further 
increased, because it supports the generation of local order. We see that the packing frac- 
tion at which the correction starts to increase strongly depends on the value of A (for the 
smallest A's it does not start to increase). Indeed, such a behavior also agrees with our 
previous consideration, because the larger A the more effective is the attractive potential 
in supporting local order, so, if the increase of the correction is due to the generation of 
local order, it must increase first for larger value of A. 
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Fig. 4. Asw as a function of the well width, T* = 1. 

In summary, we have calculated the glassiness in simple liquids described either by 
a hard sphere interaction or by a hard sphere to which a weak square well potential 
has been added. For the hard sphere system we find in the whole range of the packing 
fraction between and close packing that the glassiness curves, obtained from three 
different equations of state for the hard sphere system, are practically identical. Clearly, 
we do not expect our results to be valid in all that range because of the simple fact 
that the PY and CS equations of state are only approximate and can not expected to 
hold beyond -q = 0.5. The effect of a weak attractive potential depends trivially on the 
dimcnsionless strength of the potential (the reduced temperature). Its dependence on 
the packing fraction r\ and the ratio A between the square well and the hard sphere 
radii is more interesting. Essentially, if both -q and A are small the glassiness is reduced. 
The reason is that the attractive potential just counteracts the hard sphere. When both 
parameters are large the glassiness is enhanced because, the attractive potential supports 
generation of short range order in the system. When one of the parameters is small and 
the other is large, the resulting glassiness depends on detail but our general argument 
which is supported by the figures we present, suggests that fixing one of the parameters 
and increasing the other (the increase of r\ is limited, of course) will eventually lead to 
an increase of the glassiness. 
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